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Rocket Propulsion: Lunar Touchdown 
(A Simulation) 


1 Introduction 


For this application note, a hypothetical space mission is designed to deliver a 
delicate seismic sensor to the surface of the moon. To accomplish this, a 
probe, consisting of a two-stage rocket (a parachute will not work!) is used to 
ensure that the sensor at touchdown will arrive unscathed, with a speed of 0 
m/s. During Phase I, the probe descends from an initial orbital speed to a 
speed of -3 m/s as it reaches an altitude of 16 meters. At that point, the 
second-stage rocket is released with its engine activated for the final 7 
seconds of descent (Phase II). Just as the probe touches down on the lunar 
surface, the engine must be disengaged to prevent lift off. This simulation will 
focus only on Phase II of the probe’s final 7 seconds of descent. 


For comparison, analytical, numerical method, and analog computation 
solutions of the probe’s speed will be determined and displayed in table 1. 


2 Mathematical modeling 
Starting with the basic rocket equation (derived in Application Note #9), 
jim (Fext) = mdv/dt + udm/dt, where 
jim (Fext) = -weight of rocket, including fuel = -mg. Thus, 
-mg = mdv/dt + udm/dt 
dv/dt = -g- u/m dm/dt 
Letting dm/dt = constant = -c, where = c > 0 and m(0) = mo, 
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m= -ct + Mo 
Assuming mo = 0.900 kg, c = 0.075 kg/s (exhaust ejection rate), u = 17.2 m/s 
(engine exhaust speed relative to engine), g = 1.625 m/s? (lunar surface 
acceleration), and omitting units for simplicity, 
dv/dt = -1.625 + 1.290/(-0.075t + 0.900) with v(0) =-3.000 m/s_ (1) 
Equation (1) is used for computer setup. 


Now, for an analytical solution (fractions used for convenience): 


Multiplying the numerator and denominator of the fraction on the right by 
-40/3 and replacing 1.625 with 13/8, 


dv /dt = -13/8 — (86/5)1/(t— 12) 


dv = -13/8 dt — (86/5) dt/(t — 12) 


Vv t t 
| w= -13/8 | at—86/5 | aty(t- 12) (2) 
-3 0 0 


Integrating (2) by inspection, 
v=-3-13/8t- 86/5 In(1 -t/12) (3) 
For this simulation, the computer setup is not designed to determine the 
probe’s distance from the point of ignition to the lunar surface, with 
s(0) = 0. Nonetheless, an analytical solution will be developed. 
Letting v = ds/dt, 
ds/dt = -3 -13/8 t- 86/5 In(1 - t/12) 


ds = -3dt — 13/8 tdt — 86/5 In(1 - t/12) dt 


Ss t 


| as - sf ae—13/8 feat —26/5 fn(1-)a (4) 
0 0 


0 0 


Integration details for evaluating 
| in(1-=) at 
cae 
Letting p = (1 —t/12), then dt = -12dp. 


[im (1-=) dt = -12 J in (p)dp 


Letting u = In(p) and dv = dp, then du = dp/p and v= p. 


Integrating by parts, 
| wav = w= | vau 
| In (p)dp = In (p)p — | pdp/p 
| In (p)dp = pln(p) — | dp 
| im @yap = pant) - 1) 
Yielding, 


fro(1-L)ae=—aarfin(x-)—2) 4 


Integrating (4) in part by inspection and using the above result, 


s = -3t— 13/16 t? - 86/5(-12)[(1 - t/12)(n(1-t/12) -1) + 1] 
s = -3t— 13/16 t?+ 1032/5[(1 -t/12)(In(1-t/12) -1) + 1] (5) 
Setting v = 0 in (3) and solving for touchdown time (t’), 
0 = -3 -13/8t’ — 86/5 In(1 - t’/12) 
Using graphing techniques and inserting units, t’ = 6.594 = 7 seconds. 
Substituting 6.594 into (5) and inserting units, s = -15.838 = -16 meters. 


3 Computer setup 
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Figure 1: Computer setup for Rocket Propulsion: Lunar Touchdown 
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4 Numerical Method 
(Modified Euler method using a hand-held programmable calculator) 


PROGRAM:TOUCHDWN 
:ClrHome:ClrDraw 
"DV/DT= -1.625+1.290/(-0.075t+0.900)” 
































“INITIAL CONDITION:” 

:(0>T:-3>V 

:”STEP SIZE:” 

:0.100—-H 

Fix 2 

:Lb] 1 

‘If T>7:Then 

:Goto 2:Else 

:Disp {T,V} 

:-1.6254+1.290/(-0.075t+0.900)—>F 

:‘T+HT 

:-1.6254+1.290/(-0.075t+0.900)—S 

((F+S)/2A 

:V+AH>V 

‘Pause 

:Gotol 

:Lb] 2 

‘End 

5 Results 
Oscilloscope Time Analog Computation | Analytical | Numerical Method 
t(s) _ engine status v(m/s) v(m/s) v(m/s) 
0.0 engaged -3.0 -3.00 -3.00 
1.0 engaged -3.0 -3.13 -3.13 
2.0 engaged -3.0 -3.11 -3.11 
3.0 engaged =29 “2.93 22,98 
4.0 engaged -2.5 “2,03 22.08 
5.0 engaged -1.8 -1.85 -1.85 
6.0 engaged -0.9 -0.83 -0.83 
6.6 disengaged 0.0 0.00 0.01 
7.0 if not disengaged | 0.6 0.68 0.68 

















Table 1: Solution Comparisons 
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Figure 2: Probe’s speed as a function of time (- indicates downward) 


*For this application note, the display was produced during a single run by a 
differential equation analog computer prototype using discrete components 
with tolerances between 1% and 10%. 
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